Rapid growth of cloud droplets by turbulence 



V. Dallas and J. 

Department of Aeronautics, Imperial 



c. 



Vassilicos 
London, 



SW7 2AZ, UK 



O 

(N 

> 
O 



>> 
i 

q: 

o ■ 

CZ5 ■ 
>Y 

^ : 

Ph. 



> 
00 

in 
o 



Abstract: Assuming perfect collision efficiency, we demonstrate that turbulence can initiate and 
sustain rapid growth of very small water droplets in air even when these droplets are too small 
to cluster, and even without having to take gravity and small-scale intermittency into account. 
This is because the range of local Stokes numbers of identical droplets in the turbulent flow field 
is broad enough even when small-scale intermittency is neglected. This demonstration is given for 
turbulence which is one order of magnitude less intense than typically in warm clouds but with a 
volume fraction which, even though small, is nevertheless large enough for an estimated a priori 
frequency of collisions to be ten times larger than in warm clouds. However, the time of growth in 
these conditions turns out to be one order of magnitude smaller than in warm clouds. 



The rapid growth of warm (ice-free) cloud droplets 
from 15fim to about 50/^m in a short time, typically half 
an hour, is a well-known phenomenon which still defies 
explanation. This phenomenon is essential for under- 
standing rain initiation and cloud albedo 0, Q • 

Condensation dominates cloud droplet growth till 
about 10 to 15 fim and may be producing a narrow 
droplet size distribution [l], y|, H|- If so, a subsequent 
growth mechanism involving similar-sized particles is re- 
quired to make the droplet size distribution wide enough 
for different-sized gravitational collisions to be effective 
in a final droplet growth mechanism. Such a gravita- 
tional settling mechanism may dominate from 30 to 50 
micrometres and above 0]. The identification of what 
triggers the intermediate growth between I5(im to about 
50/im and the rate with which it proceeds are major chal- 
lenges of cloud physics. 

Various authors have already proposed that turbulence 
may be the basis of the intermediate mechanism whereby 
initially monodispersed droplets can grow via turbulence- 
generated collisions 0, H|. This mechanism is the spe- 
cific concern of the present paper. It has been suggested 
that turbulence in clouds may be generating preferen- 
tial concentrations (clustering) of droplets which would 
cause a sharp increase in collision and coalescence events 
and therefore a fast growth of droplet sizes [H, @, 0| • It 
has also been suggested that caustics may activate such 
fast droplet growth However, measurements sug- 
gest Q that, in warm clouds, the droplet response time 
t p is much too small compared to the smallest (Kol- 
mogorov) time scale r n for any significant preferential 
concentration or caustics to be observable and meaning- 
ful. Indeed, Saffman & Turner [l(| considered the case 
where the Stokes number t p /t^ is effectively zero but 
the droplet size is finite and much smaller than the Kol- 
mogorov micro-length-scalc r/ (the smallest length scale of 
the turbulence). They showed that, in this case, droplet 
sizes do not grow fast enough to explain cloud dynamics 
and statistics. Their assumption on the droplet size is 
accurate as r\ ~ 1mm and droplet radius a p ~ 10~ 2 mm 
in clouds before rain initiation and it is reasonable to 
assume, as they did, that such very small droplets are 



spherical. However, they also assumed that turbulent 
velocity gradients are statistically gaussian, and this is 
known not to be true. Small-scale turbulence is intermit- 
tent and the turbulent velocity gradients are increasingly 
non-gaussian as the Reynolds number increases 
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The neglect of intermittency may be a significant 
shortcoming because rare but powerful turbulence events 
could cause neighbouring droplets to collide and coalesce 
faster than one would expect from a consideration of the 
average properties of the turbulence field. These colli- 
sions could generate a few large droplets with high mo- 
mentum which could trigger a chain of successive colli- 
sions when travelling and falling fast through the field of 
much smaller droplets. In principle, such a chain reaction 
could lead to significant droplet size growth. Kostinski 
& Shaw 12| have already argued that rare but powerful 
events are required for droplet growth, and that these 
events may have their cause in small-scale turbulence in- 
termittency. Ghosh et al. [HI have argued that such 
rare but powerful small-scale eddies can also selectively 
increase settling velocities and thereby further enhance 
droplet size growth rates. 

In this paper we show that turbulence can generate 
fast droplet size growth without the need for small-scale 
intermittency and differential gravitational settling, even 
when droplets are too small to cluster. We place our- 
selves in a situation close to but different from Saffman 
& Turner's Close in the sense that we assume gaus- 
sian statistics of turbulent velocity gradients, spherical 
droplets of finite size much smaller than rj and a particle 
response time t p which, as a result, is very small. How- 
ever our analysis differs from that of Saffman & Turner 
flo| in that it takes into account the broad range of local 
micro-time scales of the turbulence and therefore a broad 
range of local Stokes numbers in the flow. Nevertheless, 
all these Stokes numbers are predominantly too small for 
any reasonable level of clustering to be clearly present. 

We follow the approach of Reade & Collins |l4[ who 
used a Direct Numerical Simulation (DNS) of the incom- 
pressible Navier-Stokes equations to generate the turbu- 
lent velocity field in which they integrated trajectories 
of very small spherical particles with high mass density. 
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They then applied the algorithm of Sundaram & Collins 
fl5| to simulate droplet collisions and coalescence. How- 
ever, their DNS was of 3D homogeneous isotropic tur- 
bulence without well-defined inertial range (Re\ = 55) 
and the initial Stokes numbers of their droplets ranged 
between 0.2 and 0.7 which is large enough for clustering 
to occur (see Fig. 1 in [l6| and Fig. 1 in |17j|). In- 
stead, our Stokes numbers are initially all well below 0.1 
and we opt for a DNS of 2D inverse-cascading turbulence 
which ensures a wide —5/3 energy spectrum, much wider 
than can be achieved by 3D DNS. Also, it is well known 
that the velocity gradient statistics of 2D homo gen eous 
isotropic turbulence are approximately gaussian [18l |. In 
other words, there is no small-scale intermittency in the 
turbulence we are using. 

We endeavour to use a set of conditions as close to 
real warm clouds as possible within the extent allowed 
by our computational capabilities. The ratio of the outer 
to the inner length-scales of our 2D turbulence is at least 
L/rj = 110.8. This is way below the ratio 10 5 in warm 
clouds, but high enough to have a wide inertial range (see 
Fig. [1]). The rms turbulence velocity u' is set at a value 
which ensures that if our numerical values of 77 and 
are taken to correspond to 1 mm and 10 _1 — 10~ 2 s re- 
spectively as in warm clouds, then u' is smaller than the 
usual value in warm clouds which hovers around Im/s. 
In fact, our value of u' is an order of magnitude smaller 
which makes it harder for our simulations to produce col- 
lisions and also allows the adoption of small Stokes num- 
bers, much smaller than 0.1 as in warm clouds, without 
having to take inordinately small time steps in our time- 
integrations. 

In warm clouds, the ratio of the droplct's water density 
to the density of the surrounding air is p P / Pf ~ 10 3 and 
the ratio of the mean droplet radius to the Kolmogorov 
micro- length scale before rain initiation is (a p ) jf] ~ 10~ 2 . 
In our simulations we take p p /pf = 10 3 and initial 
{a p )/rj = 0.0135 where 27r/?7 is taken to be a large multi- 
ple of the filter wavenumber k c (see Fig. [J). The larger 
we chose this multiple to be, the smaller (a p ) becomes 
and the larger the number of droplet trajectories needed 
to be integrated if we want to keep a realistic droplet 
volume/area fraction. The droplet volume fraction is 
4>v ~ 10 -6 in warm clouds but we need to set an equiv- 
alent area fraction (j>A in our 2D simulations. We do 
this by requiring that the geometrical probability of in- 
terception along straight lines between two droplets is 



the same in 2D and 3D. Given a number 7V2 of homoge- 
neously distributed droplets of radius a p in an area L^ ox 
in 2D space, a rough estimate of this 2D probability is 
4 ■ 3 ■ 2a p /(2irl 2 ) = I2a p /{irh) where l 2 = (i 2 OT /^) 1/2 is 
the average distance between droplets (about 4 droplets 
at a distance I2 from each droplet). Given a number 
N3 of homogeneously distributed droplets of radius a p 
in a volume L\ ox in 3D space, a rough estimate of this 
3D probability is 67r(2a p ) 2 /(47r£ 2 ) = 6(a v /l 3 ) 2 where 
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FIG. 1: (Color online) Energy spectrum compensated with 
A; 5 / 3 and obtained from a 2048 2 inverse-cascading 2D turbu- 
lence using the numerical method and high wavenumber filter- 
ing of [IsJ. The low-pass filtering operation and its wavenum- 
ber k c are clearly visible on the right side of the plot. This 
operation removes the scales directly affected by the small- 
scale forcing needed to give rise to a stationary turbulence 
with a A; -5 / 3 energy spectrum at lower wavenumbers k. 

it is I3 = {Llox/Nz) 1 / 3 which is now the average dis- 
tance between droplets (about 6 droplets at a distance £3 
from each droplet) . Equating the two probabilities yields 

cf>A = j (ff) 2 ^ 3 3 ' Equivalently, this means that the 
number of particles which we need to immitate a certain 
3D volume fraction 4>y in our 2D simulations is 

Our parameters are tabulated and compared below with 
the parameters of typical clouds. 



The initial particle size distribution in our simula- 
tions is a narrow log-normal with a mean particle ra- 
dius (a p ) = 0.0135?7 and a width Sa p between smallest 



and largest size such that (a p )/5a p is about 10 or larger 
(we tried up to 50 and did not find any differences in 
our conclusions). The Kolmogorov time scale f n is an 
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N Ljr\_ u_ Tr, L/(uT n ) St { a p)/v 4>v <Pa Pv/pf 

Clouds - W Urljs OTs W 10~ 4 - 10"* W" 2 - 10 3- 

DNS 2048 2 110.8 1.24 0.0074 26.7 0.04 0.0135 - 10~ 4 10 3 



TABLE I: Comparison of parameters in typical warm clouds and in our DNS. The particle response time r p = t^(-^-) 2t v> an d 

we define the Stokes number St — ^j(-^p) 2 5 where the brackets are now an average over all particles/droplets. (In keeping 
with p p and p/ which are mass densities in a volume, these expressions for r p and St are for 3D spheres, not 2D disks. It is 
more important to keep a realistic dependence on a p than 2D consistency in our model.) 



average time scale determined by the average turbulent 
kinetic energy dissipation rate per unit mass. This aver- 
age dissipation rate is proportional to 2{tr(s 2 )) where s 
is the strain rate tensor and the brackets are an average 
over all space. Hence, r„ = 1/ v / 2(tr(s 2 )). The local 
micro-time scales t v of the turbulence are determined in 
the same way but in terms of the local turbulent kinetic 
energy dissipation rates per unit mass which are propor- 
tional to the local 2tr(s 2 ). In other words the relevant 
local micro-time scales are determined by the local strain 
rates such that t v = 1/ \j2tr{s 2 ). The Probability Den- 
sity Function (PDF) of all Stokes numbers t p /t v in our 
flow is given in Fig. [2 The red solid, blue dashed-dotted 
and black dashed lines on this plot mark the values of 
t p /t v which equal {t p )/(t v }, {t p /t ti ) and (t p ) /t v respec- 
tively (the brackets being averages over all particles or all 
space accordingly). We checked that the PDF of t p /t v is 
very similar to the PDF of 1/ r r; as expected from the fact 
that the PDF of t p is very narrowly peaked, and we also 
checked that the PDF of tr{s 2 ) is peaked at zero, which 
agrees with the observation that the PDFs of t p /t v and 
of 1/t v both vanish at zero. Finally, we verified that the 
PDFs of partial derivatives of velocity components with 
respect to spatial coordinates are approximate gaussians 
peaked at 0. 

Following Reade & Collins [IH, droplets are modelled 
as very small but dense circular inertial particles sub- 
jected to Stokes drag and impulsive forces resulting 
from collisions between ith and the jth particles. These 
particles are in fact point particles in the simulations 
evolving according to 




FIG. 2: (Color online) PDFs of r p /r v at t = and at the 
times of the first 10 collisions for the case where <f>A = 10~ 4 
and rj/(2-K /k c ) = 1/5, see Tabled] The red solid, blue dashed- 
dotted and black dashed lines indicate where {t p )/(t v ) = 
2.13 • 10~ 2 , <t p /t„) = 3.67 ■ 10~ 2 and {t p )/t v = 4.03 ■ 10~ 2 
are, respectively, on the plot. The pairs of red crosses and 
green plus signs on the plot correspond to early collisions and 
indicate, on the abscissa, the values of the local Stokes num- 
bers t p /t v of the two colliding particles one time step before 
collision. The ordinate value of these crosses and plus signs is 
arbitrary except for the order of collisions, the red cross and 
the green plus sign corresponding to the first (tenth) colli- 
sion having the lowest (highest) ordinate value. The fact that 
crosses and plus signs for the same collision are always very 
close to each other reflects the fact that the first ten collisions 
are for similar-size particles. 
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the positions, velocities, masses and response times of 
particles of radius a l p , and u(x l p ,t) are the carrier fluid 
velocities at x p (t) at time t. 

A method based on molecular-dynamic-simulation 
strategies [ljl [2(| checks for collisions between particles. 
Collisions are enacted in a random order whenever two 



point particles get so close that 



\x l p -x p \\ < (a p + a p ). 



(4) 



It is commonly accepted that the coalescence efficiency 
is close to one for droplets with radius less than lOO^ni 
because their ratios of inertial force to surface tension 
are small 2l|. However, collision efficiencies may be es- 



timated to be of the order of a few tens of percent for 
droplets of 15/mr radius or so [22|, [23| even if we take 
into account turbulent enhancement of these efficiencies 



2 lj | . Nevertheless, as an initial simplifying assumption, 
we take collision efficiencies to be 1 in this study. All our 
collisions are therefore assumed to give rise to a coagula- 
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tion event thus resulting in a new circular particle with 
the following properties 



"p — '"p ' "*p 

(rripVp + m p vp)/ni}p r 



(5) 
(6) 



The positions and velocities of the remaining (non- 
colliding) particles are then advanced using 4/j/i-ordcr 
Runge-Kutta algorithm. To ensure accurate integrations 
of particle trajectories, the timestep At used in our DNS 
satisfies At <C min^r*) <C rj/u' < t v < (r n ). This 
sets a lower bound to our choice of rj and therefore a 
lower bound to the volume fraction <f>y which we can 
emulate via equation (1) for a certain number N2 of ini- 

Li23L L 2n/k c _n_ and 



tial droplets. Indeed, 



L 2Tv/k c rj (o p ) 

in our simulations ^ff- - 

<o p ) 



25.65, 



22.18 



0.245 •"«•««) 2ir/k a 

and -t^t ~ 74.1 for the initial distribution of droplet 

(Op) v 

sizes. For equation (1) to give us a realistic volume frac- 
tion 4>v ~ 10 -6 would require either extremely small 
initial N2 or extremely large 2?r ^ fcc which our time step- 
ping does not allow. (A very large 2 ^ fcc would also re- 
quire a very much larger DNS with a very much wider 
range of excited scales to be justified.) Furthermore, very 
small initial values of N2 require much longer integration 
times than we can afford. We have therefore opted for 
the following three sets of simulations: N2 ~ 1.5 • 10 6 
and r 1 /{2w/k c ) = 1/5 (case of Table I); N 2 « 5.5 • 10 5 and 
77/(27r/fc c ) = 1/10; N 2 ~ 1.25- 10 5 and r]/{2%/k c ) = 1/15. 
In the first case, equation (1) gives 4>v ~ 10~ 3 , and in the 
second and third cases equation (1) gives <j)v ~ 2 • 10~ 4 
and 4>v ~ 3 • 10~ 5 respectively. Hence, we always over- 
estimate typical volume fractions in warm clouds by two 
to three orders of magnitude, but we also underestimate 
the turbulence intensity by a factor of 10. Defining an 
a priori frequency of collisions as u'/Lb ox times the ge- 
ometrical probability of interception introduced in the 
text leading to equation (1), this a priori frequency is 
fo = — }2±v_ m our simulations and, for the case corre- 
sponding to <f>v = 10~ 3 , turns out to be 10 times smaller 
than in 3D warm cloud conditions where the volume frac- 
tion is 4>v = 10~ 6 and the turbulence intensity is 10 times 
larger than here, i.e. « WfiF . 

In the case where 77/ (27r/fc c ) = 1/5 and <py = 10~ 3 , we 
integrate the trajectories of about 1.5 million initial parti- 
cles for 22 outer time-scales tl = L/u' of our flow. Based 
on the parameters of Table U which are for this case, and 
assuming that 77 is about 1 mm and t v is about 0.1 s, this 
total integration time corresponds to less than a minute, 
which is extremely short compared to the typical fifteen 
minutes to half hour usually required for droplets to grow 
from about 15/mi to about 50/um in warm clouds. Note 
that 22Ti/2 is of the same order as fifteen minutes mul- 
tiplied by ff c . 

Numbers of collisions as a function of time arc plot- 
ted in Fig. [3J We consider two different types of initial 



conditions. One where the point particles are randomly 
distributed, in which case collisions occur immediately 
after t = because a sizeable number of point particles 
find themselves, initially, close enough for (4) to hold. 
And one where the point particles are distributed on a 
regular lattice so that (4) does not hold initially for all 
pairs of particles. It might be interesting to note (see 
Fig. [3]) that, whilst the initial evolution of the number 
of collisions is very different in these two cases, they con- 
verge towards a similar time dependence at large enough 
times t/ri = tu' / L. 

In the case where particles are uniformly distributed 
at t = 0, the first collision does not occur immediately, 
but at a time equal to about 5. It,, sa 2.67(r r? ). This col- 
lision is between particles that have a local Stokes num- 
ber T p /r l; which is between 2 and 4 times larger than 
the average Stokes number, depending on the way one 
choses to estimate it, see Fig. [2] All ten first colli- 
sions occur between time t = 5.17,, « 2.67(r^) and time 
t = 8.8t v « 4.61 (t v ) after t — and involve pairs of 
particles with Stokes numbers well above any estimate of 
the average Stokes number. They therefore occur within 
about a second and are caused by extreme events within 
the air turbulence where, locally, the Stokes number is 
higher than average. Considerations based only on the 
average Stokes number would miss these collisions and 
would therefore also miss the process initiating droplet 
size growth. Note, in particular, that the local Stokes 
number characterising the first and fifth collisions are re- 
spectively 0.08 and 0.1, just about large enough for mea- 
gre signs of clustering to be present if the average Stokes 
number had such a value (see Fig. 1 in [l(| and Fig. 1 
in0). 
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FIG. 3: (Color online) Numbers of collisions of pairs of parti- 
cles as functions oH/tl = tu /L for the case where <f>A = 1CP 4 
(see Table [TJ> . (a) Evolution up to f/n = 11 for two dif- 
ferent initial distributions of particles, (b) Evolution up to 
t/TL ~ 22 for the initially uniform distribution of particles. 

Figure|3]shows PDFs of particle radii at different times. 
At t = 0, this PDF is very sharply peaked around 
O.OI3577. By t = 3tl a second peak has appeared at 
a value about 1.4 times larger than 0.013577, and as time 
progresses more peaks appear at progressively higher val- 
ues of the particle radius. By the end of the simula- 
tion, i.e. when t/TL = 22, the largest particle has a 
radius nearly three times larger than the initial particle 
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radii. This result of our simulations is noteworthy be- 
cause </t£ = 22 corresponds to about less than a minute 
and the factor 3 is not too far from the ratio of 50/xm to 
15/im. 
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Our results support Kostinski's & Shaw's [12j sugges- 
tion that powerful rare events can cause initiation of fast 
droplet growth by turbulence. However, our results also 
suggest that this phenomenon may not require small- 
scale intermittency and/or differential gravitational set- 
tling if collision efficiency is assumed perfect. Never- 
theless, one can surely expect selectively enhanced set- 
tling velocities [HI and/or small-scale intermittency to 
increase the number and/or power of events which can 
accelerate average growth rate and thereby perhaps out- 
weigh the adverse effect of collision efficiencies which are 
typically one order of magnitude smaller than assumed 
in this work. 

The authors are grateful to Ryo Onishi, Wojciech 
Grabowski and Bcrnhard Mchlig for reading the first ver- 
sion of this manuscript and making very helpful sugges- 
tions and comments. 
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FIG. 4: (Color online) PDFs of the particle radius a v nor- 
malised by the initial average particle radius (a p ) at different 
times; in plot (a) time 3tx, (b) time 11.6tl, (c) time 22tl. In 
plot (d), six PDFs are plotted together at six different times: 
blue D> O.Otl, green A 3. Oil, red x 8.5tl, cyan □ 11.6tl, 
magenta * YT.Otl, yellow o 22t£. 



In conclusion, our calculations suggest that turbulence 
in air which is carrying water droplets so small that their 
average Stoker number is of order 10 -2 and therefore too 
small to significantly cluster can nevertheless initiate a 
process of droplet growth via collisions and coalescence 
which is fast enough for droplets to grow from 15/zm to 
50/im within about a minute. The reason behind this 
fast growth is in the wide spread of local Stokes num- 
bers in a turbulent flow. A few local flow events exist 
in the turbulent field where the local Stokes number is 
much higher than the average Stokes number and high 
enough to cause a few pairs of droplets to collide and 
coalesce quite quickly even if small-scale intermittency 
and differential gravitational settling are not taken into 
account. However, this initiation mechanism may not 
be effective if the volume fraction is too low. When we 
took the simulated volume fraction 4>v to be smaller than 
10" 3 , i.e. 4> v « 2 • 1CT 4 and cj> v w 3 • 10~ 5 , the first 
ten collisions returned by our simulations were at sig- 
nificantly later times and not all local Stokes numbers 
in these collisions were significantly larger than average 
(6 and 3 pairs, respectively, in the 4>v ~ 2 • 10~ 4 and 
4>v ~ 3 • 10 -5 cases). Nevertheless, we stress that we 
have found a turbulence-generated droplet growth phe- 
nomenon which takes droplets from 15/im to 50/im within 
a time which, if multiplied by our estimated a priori fre- 
quency of collisions, is comparable to the time required 
for droplets to grow from 15/im to 50/xm in warm clouds. 
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